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ABSTRACT 

We discuss the dynamics of microquasar jets in the interstellar medium, with specific focus on the effects of 
the X-ray binaries' space velocity with respect to the local Galactic standard of rest. We argue that, during late 
stages in the evolution of large scale radio nebulae around microquasars, the ram pressure of the interstellar 
medium due to the microquasar's space velocity becomes important and that microquasars with high velocities 
form the Galactic equivalent of extragalactic head-tail sources, i.e., that they leave behind trails of stripped 
radio plasma. Because of their higher space velocities, low-mass X-ray binaries are more likely to leave trails 
than high-mass X-ray binaries. We show that the volume of radio plasma released by microquasars over the 
history of the Galaxy is comparable to the disk volume and argue that a fraction of a few percent of the radio 
plasma left behind by the X-ray binary is likely mixed with the neutral phases of the ISM before the plasma 
is removed from the disk by buoyancy. Because the formation of microquasars is an unavoidable by-product 
of star formation, and because they can travel far from their birth places, their activity likely has important 
consequences for the evolution of magnetic fields in forming galaxies. We show that radio emission from the 
plasma inside the trail should be detectable at low frequencies. We suggest that LMXBs with high detected 
proper motions like XTE Jl 1 18+480 will be the best candidates for such a search. 
Subject headings: black hole physics — ISM: jets and outflows — X-rays: binaries 



1. INTRODUCTION 

The interaction of AGN jets with their environments 
has been under investigation for several decades, partly 
because it is easily observable through the morphology 
of radio lobes (e.g. iMilev I [1980) and X-ray cavities (e.g. 
iMcNamara & Nulsenl 120071) . Because the most powerful 
AGNs are essentially stationary in the centers of galaxy clus- 
ters, models of this interaction typically on ly consider station- 
ary atmospheres that jets run into (see, e.g jHeinz et alj2006bl 
for recent work on jets in dynamic atmospheres). The evolu- 
tion in this case can be separated into three distinct phases: 
(1) the early momentum driven phase, where the ram pres- 
sure of the jet is significant for the dynamical evolution and 
the source evolves into an elongated structure with narrow co- 
coons, (2) the energy driven phase, where the slowed-down jet 
plasma inflates supersonically expanding lobes, excavating a 
quasi-spherical cavit y (this phase i s well described by the self- 
similar solution bv ICastor et all Il975t iKaiser & Alexander! 
ll997HHeinz et alj 1998b . and (3) the late, sub -sonic evolution, 
when the radio lobes (the reservoirs of relativistic gas released 



by jets) are in pressure equilibrium with the environment. In 
the case of AGNs, this radio plasma i s buoyant in the host 
galaxy /galaxy cluster atmosphere (see iRevnolds et al]|2002l 
for a more detailed review). 

It is becoming increasingly clear that X-ray binaries 
(XRBs) are also producing relativistic je t s over a wide range 
in accretion rate dFender & Hen dry 2000; Fe nder & Ku ulkers 
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l200UlGalloetal . 2003; F ender et alJ2004l) . The inner regions 
of these XRB jets, where their dynamics is governed by the 
atmosphere of the compact object powerin g them, are very 
simila r to the in ner regions of AGN jets dHeinz & Sunvaevl 
120031) . However. iHeinz et all (120021) showed that some criti- 
cal differences exist between XRB jets and AGN jets in how 
they interact with the larger scale environment, well outside 
of the sphere of influence of the central compact object. One 
of the differences is that, in comparison with AGN jets, the 
ISM poses a much weaker barrier to microquasar jets. This 
is because the XRB jet thrust (in other words, the ram pres- 
sure delivered by the jet) is much larger in comparison to the 
inertial density of the ISM than it is in the case of AGN jets. 

The second fundamental difference between the two cases 
is, of course, that XRBs do not reside in the centers of dark 
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matter halos with stratified gaseous atmospheres. Instead, 
they travel through regular Galactic ISM with some space ve- 
locity fxRB, set by supernova kicks and orbital dynamics. 
The velocity dispersion of XRBs implies that a significant 
fraction of these sources are moving with large (supersonic) 
velocities through the ISM, which will have important con- 
sequences for their dynamics. VLBI parallax measurements 
have shown that a microquasars can be moving with veloci- 
ties in e xcess of 100 km s" 1 with respect to the local standard 
of rest (Mi rabel et al.ll200ll) . As a result, the ram pressure of 
the ISM will act on the radio plasma released by the source, 
sweeping back the outer layers of the radio lobe. 

In this paper, we argue that this aspect has important con- 
sequences for the interaction of some XRB jets with their en- 
vironment: Instead of inflating stationary cocoons, in many 
cases they will produce trails of radio emitting plasma (made 
up of relativistic particles and magnetic fields) as they travel 
through the Galaxy. A very similar situation is encountere d 
in pulsar bow shock nebulae (ICorde s 1996; Frail et all 19961) . 
where a relativistic wind inflates a channel in the ISM through 
which the pulsar is moving, and, of course, in exragalactic 
head-tail sources, which are the direct equivalent of jets prop- 
agating into a moving medium. 

The paper is organized as follows: In S|2] we will present a 
simple parametric model of the interaction of the jets with the 
ISM. In Sj3]we discuss the consequences for particle and mag- 
netic field input into the ISM, Sj4] discussed the observational 
signatures of this interaction, and ^presents a brief summary 
of the paper. 

2. THE INTERACTION OF MICROQUASAR JETS WITH THEIR 
ENVIRONMENT 

No detailed, global models of the dynamical interaction of 
XRB jets with the ISM exist as of yet. In the absence of 
such models, the first logical step is to explore a relatively 
simple sketch of how this interaction occurs. The predictions 
we will make are based on what we know about the interac- 
tion of AGN jets with their environm ent, some insight g ained 
from the scaling relations derived in iHeinz et ail d2002l) and 
on the esti mates of the kinetic en ergy output from these source 
derived in iFender et all (120051) and IHeinz & Grirnml (I2005L 
HG05 hereafter). 

Throughout the paper, we will assume that the dynamical 
time of the large scale evolution of the source is long com- 
pared to the variability time scale of the XRB and the jet (see 
£12.11 ). This will allow us to neglect the variable nature of the 
jet power in the long term evolution of the large scale radio 
lobe and instead use a time averaged constant value. 

Based on the estim ated power W cve ~ 10 37 ergs s _1 of the 
jet in Cyg X-l from iGallo et al.l (120051) . and loosely guided 
by the limit on the average jet power of (W) ~ few x 
10 37 ergs s" 1 fromHGOS 1 , we adopt a fiducial ensemble av- 
erage kinetic jet power per source of 

W={W) = lO^ergss" 1 • W 37 (1) 

We will also use the estimated average aggregate power 
from all XRB jets of W tot ~ 5 x lO^ergss" 1 This esti- 
mate does not include the even more uncertain contributions 
from neutron star jets or optically thin radio outbursts, both of 
which will increase the total power and thus the importance 
of the effects discussed here. 

1 The average power in HG05 was derived from a limited luminosity 
range, excluding sources below the luminosity completeness limit, inclusion 
of which would bring the average down. 




FIG. 1 . — Radio map of extra galactic head -tail source 3C83.1 (NGC 1265), 
rotated by 90 degrees <Odea & Owenl 19861) . 

We will further assume that the gas is sufficiently ionized 
or sufficiently dense to treat the magnetized plasma released 
by the jet and its environment as fluids and use hydrodynamic 
arguments (see Ej3.21 >. 

We will use a set of fiducial parameters for the ISM. Since 
the ISM itself is turbulent and inhomogeneous we cannot sim- 
ply assign a mean density. We will discuss the three putative 
main phases (cold, warm, and hot) separately whenever ap- 
propriate. The ISM pressure is likely dominated by turbulent 
pressure and contains contributions from thermal and non- 
thermal particles and magn etic fields. We will use a fiducial 
mean total ISM pressure of dCoxll2005l) 

Pism ~ 3 x 10~ 12 ergscm~ 3 • p_ii.5 (2) 

We will assume throughout the paper that the working 
surface, i.e., the place where jet and ISM interact directly, 
slowing down the jet material (taken to be a strong terminal 
shock 2 ) is efficient at slowing the jets down to sub-relativistic 
velocities and thus dissipating most of the kinetic energy into 
internal energy of the radio plasma. For a jet with length ij et 
and half opening angle a = l°ai (where a% is the opening 
angle measured in degrees), this is the case for 

/ jct > 10 16 cm • J™*Ll- (3) 

V Til Oil 

where we define the ISM particle density as n\ = 
^ism/1 cm" 3 . The size scales we are interested in are much 
larger than this, so we are justified in assuming the working 
surface dissipates most of the kinetic energy. We will also 
assume that this plasma does not mix with the ISM due to 
magnetic flux freezing (see, however, ^3 . 2b . 

2.1. Early evolution 

During the early lobe evolution the ISM pressure Pism> 
buoyancy forces, and the ISM ram pressure Pismire c ' ue 
to the space velocity i>xrb of the XRB are negligible com- 
pared to the ram pressure p\su_v\ on the lobe due to its 
radial expansion velocity ur . It is then appropriate to use 
the well-known 1 -dimensional energy-driven bubble model 
dCastor et alj 119751: iKaiser & Alexander! 119971: IHeinz et alj 
Il998t iGalloetalJ 120051) as sketched in Fig. |2 This model 
assumes that the XRB releases relativistic plasma in a spheri- 
cally symmetric fashion 3 

Assuming (a) an adiabatic equation of state, (b) uniform 
pressure for the radio plasma inside the lobe, (c) ram pressure 

2 There is little fundamental difference if instead the jet material is decel- 
erated through a series of weak shocks along the jet 

3 This assumption is motivated by the fact that the internal sound speed 
in a relativistic plasma is much larger than the expansion velocity of the ra- 
dio lobe. Thus, regions of overpressure (i.e., the working surfaces of the 
jets, where the jets interact with the ISM) will quickly distri bute their over - 
pressure within the relativistic plasma cocoon (e.g. Begelman & Cioffi 1989). 
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FIG. 2. — Sketch of early and late evolutionary phases in XRB radio lobe 
dynamics. Note that in the early stages the expansion velocity of the plasma 
is much faster than the space velocity of the XRB, thus the evolution is quasi- 
spherical. Later, the ram pressure due to the XRB space velocity and the 
thermal pressure of the ISM dominate, and the stripped plasma leaves a trail 
behind the XRB. 



balance with the ISM, and (d) energy conservation, the Castor 
solution is 

*3\ i/5 
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■early 



0.65 



(w)t 3 

PlSM 



(4) 



Apart from the expected synchrotron emission from inside 
the radio lobe, the supersonic expansion into the ISM will 
also lead to therma l emission from the shocked ISM shell 
dKaiser et al.l 120041) . similar to expanding supernova rem- 
nants. Shocks are extremely valuable as diagnostic tools, 
since the shock strength measures the expansion velocity of 
the source, which, when coupled with the source size, can tell 
us the age and average kinetic power of the source. An ex- 
ample of such a source is Cyg X-l, for wh ich this thermal 
emission has been detected (Gallo et al. 2005|)- 

2.2. Late evolution: Head-tail sources 

At sufficiently late times (the case we are interested in here) 
the jets will have traveled sufficiently far for the dynamical 
pressure due to the XRB motion to become important. 

The dynamical pressure of the XRB motion through the 
ISM with velocity vxrb will become comparable to the jet 
thrust (i.e., the momentum flux along the jet) when the ad- 
vance speed of the working surface equals the space velocity 
of the source, «head = «xrb- If this condition is satisfied, the 
working surface can no longer advance in the direction of the 
jet and will thus become stationary in the frame of the XRB. 

we 



Given our fiducial sound speed of c soun d 
will express the XRB's space velocity as 



= 17kms~ 



«XRB 

17kms _1 



(5) 



Furthermore dynamical instabilities and precession are believed to spread the 
jet thrust out over a much larger solid angle than subtended by the opening 
angle of the jet for the case of AGN jets I Scheuei 1982), driving the aspect 
ratio of lobes toward unity. We will assume that the effective half opening 
angle is 9 = y^Q/ir ~ 20° (e.g. lYoung eTaI]|2003) where H is the solid 
angle swept out by precession and jitter of the jet axis. 



The shock jump conditions and pressure balance across the 
working surface (the terminal shock) imply that this happens 
when the jet thrust 



Pjet 



5T + 3 W\ 



jet 



r Vjot i2v{eiy 



equals the total dynamic pressure 



Pdyn = 4.1 prSM^M 

TZ M = (9 + 10M 2 [l 



(6) 



(7) 



+ v 7 ! + 9/(4M 2 )] } / 



/37 



of the ISM, where M = «xrb /c SO und is the effective Mach 
number of the XRB's motion through the ISM with re- 
spect to the effective sound speed of the ISM, c soun d = 
\/5pism/3pism- We adopt a fiducial value of T = 5T 5 for 
the bulk Lorentz factor of the jet, which is, however, of little 
importance throughout the rest of the paper. 

The expression for pdyn in ec l- can easily be derived 
from the sh ock jump conditions on eac h side of the working 
surface (e.g. Blandford & McKee 1976), with non-relativistic 
equations of state for the unshocked cold jet and ISM gas 4 . 

The two dynamical pressures pd yn and pj Ct are equal when 
the length of the jets reaches 



Ijct — 



6 x 10 18 cm 



W37 25T + 15 
p_ii.5 28173ft 



M 



1/2 



(8) 



where the expression in brackets is unity for M = 1 and the 
fiducial parameters and j3 — Wj e t /c. At this length, the work- 
ing surface becomes stationary in the frame of the XRB, and 
since the ISM is moving at velocity vxrb with respect to the 
XRB, the plasma that is released at the working surface must 
be swept back away from the XRB with the ISM at velocity 
^xrb an d create a trail of "debris" plasma, left behind by the 
XRB jets. 

The time for the jet to propagate that far into the ISM is ap- 
proximately t cq w Ijct/ vxrb ~ 10 5 yrs, which i s long com- 
pared to typical accretion disk time scales (e.g. iFrank et all 
|2002|) but still short compared to typical companion (i.e., 
XRB) lifetimes. It is therefore necessary to further discuss 
this quasi stationary phase of the large scale dynamics of ra- 
dio plasma released by XRBs (it is quasi stationary since the 
XRB space velocity vxkb an d the local ISM density tiism 
change as the XRB traverses the Galaxy). 

2.3. The debris: Radio trails 

As mentioned above, we assume that the working surface 
is stationary in the frame of the XRB. Its propagation velocity 
into the ISM must therefore be i>xrb and its pressure must be 
given by Pd yn from above. 

The energy released by the jet at the working surface must 
then equal the jet power: 



W iei = 



7 dV head 

'Pdyn " 



r 



dt 



(9) 



where Vhead is the volume of plasma streaming out of the 
working surface. The factor [7 a d/(7ad — l)]Pdyn is the post- 
shock enthalpy of the jet plasma and includes both the internal 
energy pdV/(^ a d — 1) of the plasma released and the pdV 

4 The 7 = 5/3 equation of state is not a crucial assumption and can 
easily be modified to reflect a relativistic equation of state of the jet gas (for 
example). The coefficients only change slightly and the qualitative results are 
unaltered. 
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work done on the ISM. Since 7 ac i = 4/3 for the jet plasma, the 
volumetric rate at which jet plasma is released at the working 
surface is 



head 



Hi 



jet 



dt 



4pdyn 



(10) 



Since far downstream (along the swept back trail of radio 
plasma away from the source) the pressure has to equal the 
ISM pressure (otherwise the energy in the Mach cone behind 
the XRB would diverge), the radio plasma originally released 
at pressure pa yn must expand adiabatically (thus doing further 
work on the ISM). Thus, the asymptotic volumetric rate at 
which plasma enters the trail is 



dV_ 



Wjct ( Pdyn 
4pdyn VPISM 



jet 



4pdyn 



Vpism/ 



3/4 



(ID 



In order to estimate the asymptotic lateral radius R cq of the 
trail (approximated as a cylinder around the XRBs trajectory) 
from its inflation rate (eq. [TTJ, we need to know the asymp- 
totic flow velocity i? tr of radio plasma through the trail. The 
pressure gradient behind the bow shock (from the ISM ram 
pressure pdyn to pism) will accelerate the plasma along the 
channel. 

In principle, it is possible that the plasma will pick up 
most of its original outflow velocity and reach mild ly rela- 
tivistic speeds. For example, Buccianti ni et alj ((2005) inves- 
tigated the dynamics of pulsar bow shock nebulae using axi- 
symmetric 2.5D relativistic MHD simulations and found that 
the velocity inside the channel reaches values up to 0.6c. They 
concluded, on the basis of their simulations, that the channel 
width should be comparable to the standoff distance between 
the pulsar (or XRB in our case) and the stagnation point at the 
head of the shock. However, the enforced symmetry in their 
model makes the solution they found over-stable. In reality, it 
is likely that dynamical instabilities (e.g., Kelvin-Helmholtz 
or kink instabilities, depending on the strength and topology 
of the magnetic field) would disrupt such rapid flow through 
the channel. 

We expect the average flow velocity to be roughly that of 
the ISM sheet around the channel. Asymptotically, this veloc- 
ity must be UxRB (otherwise the total kinetic energy would di- 
verge). Thus, we will assume that the plasma, far downstream 
from the XRB will have a flow velocity of v tI = t>xRB V with 
V > 1. 

Note that even in the case of V ^> 1, the plasma must still 
come to rest in the ISM at some point, where it would inflate 
a bubble. This would be around the birthplace of the XRB or 
around any major change in environmental parameters (i.e., 
changes in ISM density or pressure), where the channel width 
would change significantly, leading to a series of "sausage 
link" bubbles, connected by narrow channels. Note also that 
the total volume of plasma released should not be affected. 

The equilibrium radius R cq of the roughly cylindrical trail 
is then 



Ron 



dV 



1 



dt -kvxrbV 



>2 x 10 20 cm • 



.37 



P-11.5«17 V 



7? 8 



(12) 



This is much larger than the value of lj Ct (eq. |H). Thus, 
the ultimate cross section of the trails is much larger than the 
equilibrium length of the jets and our assumption that both jets 



are feeding a single trail is justified (this is not qualitatively 
important for our argument, however). 
It is worth pointing out the implicit weak dependence of 

— 1/8 

R eq on density through 1Z M which implies that the size 
scale of the trails produced should be similar in different en- 
vironments - an increase in density by three orders of mag- 
nitude (with otherwise identical parameters, i.e., pressure and 
jet power) only reduces R cq by a factor of 3.5, while a de- 
crease in density by three orders of magnitude (i.e., a hot ISM 
phase environment) only increases R eq by a factor of 1.3. 
Thus, the most significant effects on R cq are likely to come 
from W 3 7.3, P-n.5, and v 17 . 

For a binary lifetime of t = tg x 10 8 yrs (where tg is the 
age of the source in units of 10 s yrs), the length of the trail 
left by the source is 



trail 



= 5 x 



10 21 cm 



^17*8 



(13) 



Based on this expression, a source will produce a significant 
plasma trail if I ^> R cq , or 



1/2 3/2 1/2^-1/2 | 
-11.5 U 17 v vv 37 /C M 



> 0.06 



(14) 



If this condition is not satisfied, the source dynamics are better 
approximated as spherical and stationary (if significant gradi- 
ents exist in the environment of the XRB, the source evolution 
will, of course, still deviate from a quasi-spherical, symmetric 
case, as might be the case in Cygnus X-l). 

The fiducial values used in eq. (fl4l are appropriate for 
LMXBs, implying that LMXBs typically leave trails in the 
ISM, especially the high velocity tail of the LMXB popula- 
tion with vn 3> 1. Since HMXBs have low space velocities 
and short lifetimes, this condition is typically not satisfied for 
HMXBs: For a characteristic HMXB velocity of v u V — 0.35 
and a companion lifetime of tg < 0.1, the expression on the 
left is 0.015. Thus, HMXBs typically should not produce 
trails, except for very low power source, and sources in the 
very high velocity tail of the HMXB distribution in dense, 
high pressure environments (such as the Galactic center). 

In arguments above we have implicitly assumed that a fluid 
treatment is valid, i.e., that the relativistic particles released 
by the jet interact collectively with the ISM. This is equiv- 
alent to assuming that the particles are frozen to the plasma 
and that microscopic mixing between ISM and trail plasma is 
not a dominant effect. We will discuss this assumption and its 
validity in more detail in 13. 21 Any significant amount of mix- 
ing will only make the conclusions we will draw below about 
seeding of magnetic fields into the ISM stronger (however, it 
might alter the observation appearance of trails as well). 

3. MIXING AND THE IMPACT ON THE ISM 

Given the estimates for the properties of the radio lobes and 
trails produced by the ensemble of Galactic XRBs, we can 
estimate their impact on the ISM. Before going into detail, 
we will lay out our assumptions about the ISM that the XRBs 
are interacting with: 

For the XRB velocities, we will use a Maxwellian distri- 
bution with mean velocities of oh ~ 6kms~ 1 and <tl ~ 
16kms _1 for HMXBs and LMXBs respectively, which re- 
produces the scale heights of both populations. A fraction of 
about /disk ~ 75%/disk,75 of the LMXBs reside in the Galac- 
tic disk. 

For the ISM density distributions we will use volume 
filling fractions of / hot ~ 70%, / waim ps 25%, and 
/cold ~ 5% for the hot, warm, and cold phases respectively 
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dMcKee & Ostrikerl [l977h . We will assume that LMXBs 
are randomly distributed among the three phases since they 
are old enough to have moved far enough away from their 
birth places to have forgotten about their initial environments. 
HMXBs do not live long enough and have velocities too small 
to have traveled very far from their birth place, about 60 pc 
for vxkb = 6 km/s and a life time of 10 7 yrs. This has to be 
compared to typical molecular cloud size of about lOOpc. We 
will thus assume that / co id = 100% for HMXBs to be conser- 
vative (smaller values of / co id will lead to larger values of the 
total volume of plasma released by HMXBs below). 

3.1. The total volume of radio plasma released into the 
Galaxy 

The total volume in radio trails produced by XRBs depends 
on the distribution of i>xrb an d Pism only through TZ 1 ^ 4 . 
Thus, estimates of the volume of plasma released by jets will 
be robust against uncertainties in those quantities. 

With the ISM parameters given above, the volume in trails 
produced by the LMXB disk population is simply the cur- 
rent rate of volume injection per LMXB, averaged over the 
velocity distribution of LMXBs, multiplied by the number 
of LMXBs, multiplied by the age of the LMXB population 
i L = xl0 10 ti yrs. 

Using the estimated total integrated kinetic power of the 
LMXB population of W L ps 3.5 x 10 38 W L ,38.5, this gives: 

Vlmxb .disk 



WL, 38. 5/disk, 75^10 



' 5.5 x 10 66 cm 3 • 

■ht(ch: /cold .warm, hot ) 



(15) 



where is unity for our fiducial parameters and a very 
slowly varying function of cr L , m, and / co i d ,warm,hot OVd 
the parameter range considered. The velocity dispersion <tl 
drops out of eq. ( TTBT l to lowest order. Clearly, the most un- 
certain parameter left is Wlmxb,38.5i on which Vtot depends 
linearly. 

Taking fhot = 100% for the halo population, the vol- 
ume produced by the halo LMXBs is about VL,haio ~ 2 x 
10 66 cm 3 /p_n 5. However, a sizeable fraction of the plasma 
produced in the disk will eventually leak into the halo (see 
below). Since p-n.5 is much smaller in the halo than it is in 
the disk, the total volume filled by plasma released by halo 
LMXBs and that emanating buoyantly from the disk (see be- 
low) is thus likely much larger than this estimate. 

For the HMXBs, the total volume rate should be integrated 
over the star formation rate history, assuming a current star 
formation rate of 3 Af© yr _1 M3 (where M3 is the Galactic 
star formation rate in units of 3 solar masses per year). 

With the slightly lower value of the total power Wr and 
the total star formation rate integral of J dtMspn = Meal 
ID Mq Mh (where Mn is the Galactic stellar mass in units 
of 10 11 solar masses), this gives 



t t . 1n 66 3 WH,38.3-Mil 

P-11.5M3 



M<W«0 (16) 



Again /ih is slowly varying. 

Summarizing, the total volume of radio plasma released 
into the disk by XRB jets over the lifetime of the Galaxy is 
then roughly 8 x 10 66 cm 3 , and the volume released by halo 
LMXBs is roughly 2 x 10 66 cm 2 x (p d isk/Phaio)- 



This should be compared to the disk and the halo volumes: 
Taking the disk scale height to be 500 pc and the disk radius 
to be 10 kpc, the disk volume is Vd; s k ~ 10 67 cm 3 , while the 
halo volume is roughly Vhaio ^ 10 68 cm 3 . Clearly, XRBs 
will have released sufficient relativistic plasma to fill a sig- 
nificant fraction of the disk and even the halo. If the radio 
plasma were confined to the disk, it would fill up to 100% of 
the disk (and could account for all of the hot ISM phase in 
the disk). This immediately implies that the buoyancy of this 
plasma in the vertically stratified Galactic gas will eventually 
move most of the trails out into the halo. It also implies that 
a volume fraction larger than 10% of the halo might be filled 
with radio plasma produced by XRBs. We will discuss further 
implications of this result below. 

3.2. Microscopic mixing 

Microscopic mixing of radio plasma and ISM is only go- 
ing to be important in sufficiently neutral ISM phases, since 
otherwise the strong magnetization of the trail plasma will ef- 
fectively decouple both phases. 

While the ISM will be sufficiently ionized near the working 
surface, it might well be neutral far back ablong the trail. In 
this case, neutral particles can enter the trail. The cross section 
for collisional interaction between the neutral particles and the 
cosmic rays of the trail is negligible. Thus, the only way for 
neutral ISM particles to interact at all with the trail plasma is 
through interaction with ionized ISM particles via ambipolar 
diffusion and charge exchange reactions (we will ignore the 
latter given the expected low cross sections). 

While a detailed discussion of microscopic mixing of rel- 
ativistic plasma across radio lobe boundaries is beyond the 
scope of this paper, we will briefly present an order of mag- 
nitude estimate of the importance of ambipolar diffusion of 
neutral ISM into XRB radio lobes and trails to show that the 
exptected am ount of mixin g is ast rophysically interesting. 

Following iDraine et al.l d 1 983b . we use a constant ion 
drag coefficient of 7i on = 3.5 x 10 13 cm 3 g _1 s _1 . In 
the cold neutral phase, it is appropriate to use an ioniza- 

tion fraction of & = rii/n n = Cip n with Q w 3 x 
IQ-is gi/2 cm -3/2 ( Elmegreenl 1 1 9791) . In the warm neutral 
mediu m, the ionization fraction is higher, of order £; ps 10% 
(e.g. lDalgarno & McCravll 19721) . 

Next, we need to estimate the field strength inside the trail. 
Given the estimated pressure and size of the plasma trail, we 
can parameterize the field strength in terms of its equipartition 
value, 



B cq = ^36*871-^1/(1 + 60 ( 17 > 

«H/iG • V2£bP-ii.5/(1 + £b) 
where £g is the ratio of magnetic to particle pressure and p tot 
is the total (particle plus magnetic) pressure. For simplicity, 
we assume the field is ta ngled isotropically, leadin g to the ad- 
ditional factor of 3 (e.g. lHeinz & B egelman 200(J. 

The ambipolar diffusion c oefficient f or cold neutrals enter- 
ing the trail/lobe is then (e.g. lShul[l992l) 

v\ B 2 1 

r) A 

-^amb.cold A ^ 

7Pi 47T j O n 7 6; p n 

/ 1 n4 —3 \ 3/2 

a m20 2 -1 ' 10 4 cm d x 

4 x 10 u CYTL s 1 • p-11.5 



2£e 



1 + Ce 



and for warm neutrals 

^amb,warm — 10 Cm S 



p-ii. 5 10% 2£ B 

< € i + £e 



(18) 
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We are interested in the asymtotic mixing fraction far back 
along the trail, when it has reached its equilibrium radius. 
Thus, we will not consider dynamical effects due to the trails 
initial expansion on the total mixing fraction 5 . 

A conservative assumption about the longevity of the radio 
trail is that it will last at least one sou nd crossing time against 
any dynamical instabilities (see 33.31 1, Within a sound cross- 
ing time T soun d = R e q/ Csound, the ambipolar diffusion length 
is 



^amb 



D 



amb 



(19) 



Tsound 

Thus, a lower limit to the mixing fraction of thermal material 
from the cold neutral phase into radio trails/lobes produced by 
microquasars is 



mix, cold : 



>2% 



and for warm neutrals 




(20) 



I mix, warm 



>io% 




We conclude that, for the warm and cold neutral ISM, the 
magnetic field inside the trail presents only a moderate bar- 
rier and mixing could be relatively efficient, quickly creat- 
ing a boundary layer of a few percent thickness around the 
trail/lobe. At the same time, ambipolar diffusion is probably 
not strong enough to affect the dynamics of the radio plasma 
on a global scale. 

In conclusion, mixing is significant in the neutral ISM 
phases, with mixing fractions of order of a few percent. It 
is likely orders of magnitude smaller for XRBs in ionized 
environments 6 . This suggest that the global average of f m i x 
should be about a factor of 10 smaller than the volume filling 
fraction of neutral gas in the Galaxy. So, given recent esti- 
mates of v olume filling factor of up to 40% for the warm neu- 
tral phase (Heiles & Troland 2003) we "guesstimate" a glob- 
ally averaged mixing fraction of order f m - lx < 4%. 

3.3. Competing destruction processes 

As the radio plasma expands into the ISM a number of pro- 
cesses will compete to either mix the two phases or remove 
the plasma from the Galactic disk into the halo. 

3.3.1. Buoyancy 

The ISM is vertically stratified due to the gravitational field 
of the Galaxy. The low density radio trails are buoyantly un- 
stable in this gas, which will drive them upward, away from 
the Galactic plane. 

The buoyancy time scale for a cylindrical trail of radius R 
is set by the downward ISM dynamical pressure and the up- 
ward buoyancy of the radio plasma. For a Galactic rotation 

5 During the radial expansion phase of the lobe and in the radially expand- 
ing section of the trail, ambipolar diffusion will only work as long as the 
diffusion speed is faster than the expansion velocity 

6 Anomalous diffusion can lead to significant m ixing even between two 
completely ionized phases (e.g. Narayan & Medvedev 2001). 



velocity v r ~ 220 km s _1 (which sets the vertical gravity g z ), 
a Galacto-centric distance of r ~ lOkpcrio (where rio is 
the Galacto-centric distance in units of 10 kpc), and a vertical 
trail position z ~ 200 pc Z0.2 above the Galactic mid-plane 
(where 20.2 is the trail height in units of 200 pc), and assum- 
ing a drag coefficient C\y of about 1 for a cylindrical trail, the 
buoyancy speed is 



g z TrR v rot r - 5 - 
C w r 



(21) 



;3 x 10 5 cms _1 



1/2 

Z 0.2 



no VP-11.5U17 



1^37 



1/4 



For a disk scale height of H 500 pc H , 5 , the buoyancy 
time is 



T b ~ H/ v B 



(22) 



1.6 x 10 s yrs 



rio #0.5 ( P-11.5«17 



1/2 

Z 0.2 



W37.3 



1/4 



which is, be default, longer than a sound crossing time (since 
buoyancy is always sub-sonic). This is comparable to the life- 
time of an HMXB, but could be short compared to the lifetime 
ofanLMXB. 

Unhindered, this process will eventually remove all of the 
radio plasma from the disk into the halo. However, a number 
of processes will work against this transport out of the plane, 
trying to mix the relativistic plasma into the ISM. In order 
to estimate their importance, one needs to compare the time 
scales over which they act to the buoyancy time. If they are 
occurring on more rapid time scales, they might mix a signif- 
icant amount of radio plasma macroscopically with the ISM 
before it has time to escape into the halo, especially because 
all of these processes will only increase the trails surface area 
and thus its buoyancy time. At the same time, an increase in 
surface area will also increase the microscopic mixing frac- 
tion / mix . 

3.3.2. Shredding 

Due to the differential rotation of the Galaxy, any object of 
finite size will experience Galactic shear. Since the radio trails 
are embedded in the ISM, they will be affected. 

Taking the Galactic rotation velocity v r ~ 220 km s -1 to 
be constant with radius, the angular velocity shear across the 
trail is duj/dr = —v r /r 2 . The Shredding time for an object 
much smaller than its Galacto-centric distance r that is subject 
to this shear is 



7shear 



r/vj ~5x 10 7 yrs • 



rio 



(23) 



independent of the actual size of the object. 

After several shredding times the object will have been sig- 
nificantly stretched and distorted and will eventually be mixed 
with the ISM. Note that the shredding time is approximately 
equal to the buoyancy time, so for an XRB in the Galactic 
plane, it is plausible that shredding will distort, and possibly 
even destroy the radio trail before it can escape into the halo. 

3.3.3. RT and KH instabilities 

In addition to shear-shredding, the trails will also be subject 
to dynamical instabilities. The same buoyancy process that 
drives the trails out of the Galactic plane will also trigger the 
formation of Rayleigh-Taylor instability, acting to destroy the 
trail and macroscopically mix the radio plasma with the ISM. 



Microquasars as head-tails sources 



7 



The growth time for a mode with wave number k is r ~ \fgk 
(for large density contrast between trail and ISM). Since the 
trails are filled with magnetized plasma, magnetic tension will 
suppress the growth of small wavelength modes. However, 
global modes will likely not be suppressed by tangled fields 
(at least in the linear regime), and so the growth time for a 
body mode with wave number k = 2n/R cq will be 



2ng 



R 

2ttz 



R 



1 



H V 4CV 



< r b (24) 



with Cw ~ 1 being the drag coefficient. Thus, RT instabil- 
ity will set in well before buoyancy had time to remove the 
plasma out of the disk and into the Halo. If RT instability 
does grow at the hydrodynamic rate, it will fragment the ra- 
dio trails into smaller pockets, which will be subject to further 
RT instability. 

Eq. (l24l i shows that RT instability is more efficient the 
smaller the size of the plasma pocket, thus, ;/RT instability is 
not suppressed on scales comparable to R cq by the presence 
of large scale magnetic fields, it will destroy the radio trail 
into a spectrum of bubbles small enough for magnetic tension 
to stabilize them before buoyancy transports the plasma out 
of the disk. Since the buoyancy time of smaller bubbles is 
longer, the RT fragmentation will delay buoyant transport and 
give other processes more time to act and, by increasing the 
surface area of the plasma, increasing the efficiency of diffu- 
sive processes. 

In addition to RT instability, Galactic differential rotation 
and buoyancy drift will induce KH instability along the edge 
of the trail. The typical velocities encountered are of order 

"buoy ~ v r VRz/r and u s hcar ~ 2v r R/r. Ushoar will domi- 
nate at early times, «buoy will become larger after a few buoy- 
ancy times. Taking the relativistic inertia of the radio plasma 
into account, the KH growth time of a body mode with A ~ R 
is roughly 



R 2 



bouy 



e 



Ap 



ISM 



12(^ 2 ^ hoar + z^) c sound 



(25) 



(26) 



Clearly, the ratio c/c soun d is much larger than either z/R or 
H/R for typical trail radii (see eq. fT2l . and KH instability 
should not contribute much to the global destruction of radio 
bubbles and trails. However, it could well induce a spectrum 
of short wavelength modes that creates a shear layer around 
the trail, where efficient mixing between the two phases oc- 
curs. 

3.4. The impact on the ISM: Magnetization and cosmic ray 

seeding 

Given the total volume of plasma released by XRBs and 
some value of / m i X , we can estimate the strength of the mag- 
netic field that is released and seeded into the ISM. 

We will use the expression for the field strength inside the 
trail from eq. ( TTTl ). Further assuming magnetic flux conserva- 
tion, the strength of the magnetic field after the residual frac- 
tion / m i X of the radio plasma is mixed with the ISM is 

-^mixcd ~ /mix-^eq (27) 



2£e 



h _ 5.5WL,38.5/disk,75*10ftL + 4W H ,38.3M 1 ife H /M3 

9.5 

where b is unity for the set of fiducial parameters. Since tio 
and Mn enter linearly into this expression, the magnetic field 
provided by XRB jets will have been linearly increasing as 
a function of time for most of the Galactic history, but will 
have been increasing more rapidly during intense star forma- 
tion episodes, when M increased dramatically, since then the 
HMXB population would have been much more numerous. 

This should be compared to an estimated mea n Galactic 
field of order B G ~ 10 (e.g. IWentzell fT%3h . which is 
believed to have been created by dynamo action on a small 
seed magnetic field. Comparing both values implies that 

/mixV 2 Wb-ii.5(l + £b)]& < 1. which is easily satisfied 
since we estimated f m m 0.04 above. 

The seed field required to produce the observed G alactic 
field is rather small (B sccnl ~ 1 ~ 1 7 G [Anderson! 1 9921) . Even 
under the most conservative assumptions possible — mixing 
of only one skin depth of the thickness of one proton gyro 
radius (/ m j x ~ 10 -9 ) — such a field strength could have 
been provided by XRB jets within about 10 5 yrs after the first 
XRBs turned on (though it would need a few Galactic revo- 
lutions of duration r rev ~ 2.5 x 10 s yrs to get spread over 
a significant fraction of the disk). More reasonable choices 
of /mix imply that the strength of the seed field provided by 
XRB jets is likely orders of magnitude larger. In fact, given 
the numbers derived above, is likely that the strength of the 
field seeded by XRB jets is only about two orders of magni- 
tude below the observed mean Galactic value. 

Of particular importance in this case is the fact that LMXBs 
travel far from their birthplaces, leaving behind trails of mag- 
netized plasma, which should permeate the Galactic plane 
more or less uniformly (modulated by the average space den- 
sity as a function of Galacto-centric distance). 

Since XRB production is an unavoidable part of star for- 
mation (Grimm et al. 2002|) an d since jet ac tivity is an inte- 
gral part of XRB activity (e.g. lFenderll2001l) . it is clear that 
jets must play a role in providing magnetic fields in form- 
ing galaxies and the cumulative effect of XRBs on the ISM 
implies that they are likely responsible for a non-negligible 
fraction of the field in mature galaxies as well. 

Along with the magnetization, it is clear that relativistic 
particles from within the lobe will mix into the interstellar 
medium as cosmic rays. T his process was discussed in more 
detail in lHeinz et al.l (120021) . A critical question raised in that 
paper was the adiabatic cooling of these particles. The pres- 
sure inside the working surface has already been estimated 
in eqs. (|7]i and the pressure inside the trail/lobe to which the 
cosmic ray plasma has to expand is simply the ISM pressure. 

Taking the bulk of the particles leaving the working sur- 
face to have energies comparable to the randomized specific 
kinetic energy of the jet, i.e., 7 ~ T, the final energy of any 
cosmic ray protons in the trails is 



E, 



OR W 7 TO pC 



PlSM 
Pdyn 



3.3 GeV ■ T 5 



(28) 



P-u.b(1 + 6j) 



indicating that adiabatic losses are likely small enough to 
leave the results regarding the total e nergetic of cosmic ray 
injection by microquasars discussed in lHeinz et alj d2002l) un- 
changed: the total energetics of microquasars is likely to con- 
tribute of order a few percent to the Galactic cosmic ray pop- 
ulation. Thus, microquasars located in the vicinity of molecu- 
lar clouds should be observable through gamma-ray emission 
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due to hadronic interaction with the cloud protons detectable 
by GLAST. A corollary of this process is GeV neutrino emis- 
sion that should also arise. The high energy tail from this 
emission might be detectable with high energy neutrino de- 
tectors like Ice Cube. 

4. OBSERVABILITY 

Given the estimated pressure and size of the plasma trails 
and working surfaces (i.e., hot spots), we can estimate their 
synchrotron luminosity and surface brightness. We once again 
use the expression from eq. ( fTvT > to parameterize the field 
strength. Since it is unclear what the composition of the 
plasma is (e + /e~ vs. p + /e~), we parameterize the proton 
pressure as p p = £ p p tot /[(l + £b)(1 + £ p )], where £ p is the 
ratio of proton to lepton pressure. 

4.1. Radio trails 

While the field strength we estimate for the trails is not 
much larger than the mean Galactic field, it is safe to assume 
that the partial pressure from relativistic electrons inside the 
trails is larger than outside, implying that the trails will pro- 
duce synchrotron emission that is significantly enhanced over 
the mean Galactic emission. 

For a powerlaw distribution of electrons with index 2 such 
th at dN/d'y oc 7~ 2 , the synch rotron emissivity of the plasma 
is dRybicki & Lightmanll979h 



5 x 10 d8 ergs / 2p_n 



cm 3 s Hz 



1+6 



1 + £ P V 5 GH 



(29) 



The magnetic field strength inside the trail is approximately 
B cq ps 1.4 x 10~ 5 mGp_ii. 5 26/(1 + £b), which gives a 
cooling time of 



7*cool 



10 7 yrs • 



1 + 6 



11. E 



3/2 



-1/2 



(30) 



for an electron emitting at frequency v = 5 GHz v§ . The 
length of radio trail a particle can travel before cooling losses 
affect it is 



cool 



ol^XRB 



5 x 



10 20 cm 



1 + 6 



11.5 



3/2 



1/2 



(31) 

which is comparable to R cq from eq. (fLZb . Thus, the trails 
left behind by XRBs should typically only emit at frequen- 
cies well below 5GHz. Above such frequencies, the radio 
emission from XRBs should be concentrated around the XRB 
itself. Given their large angular sizes and low frequency emis- 
sion, these trails will be ideal targets for the new wave of 
high fidelity, low frequency radio arrays, such as LOFAR, the 
LWA, the MWA, and, ultimately, the SKA-low. 

The central surface brightness of a given trail within a cool- 
ing distance Z coo i from the XRB will be 



s 2R cq e 
arcsec 2 



(32) 



1+6 



1 + 6 V^itV M 



and the integrated trail luminosity at a fixed frequency is 



3 x lO^ergsHz^s" 1 



(33) 



While every microquasar should have some kind of radio 
lobe/trail in its vicinity, some sources will be better suited to 



study this effect than others. Cygnus X-l, for example, is 
already known to possess a dark radio lobe. As an HMXB 
with a relatively young age, it is unlikely to show much of a 
radio trail, as argued above, and indeed, the radio lobe seems 
to be close to the source, though it does show a one-sided 
asymmetry: the radio ring surrounding it that led to the claim 
of the radio lobe is one-sided. This asymmetry could be due 
to the fact the source is located at the border of a high density 
cloud. 

Arguably one of the best sources to look for a radio 
trail is XTE Jl 118+480, which is known to have both a 
steady low/hard stat e jet as well as a large space velocity 
(Mira befet al.l l2001). We suggest that this should be one of 
the first targets to study with low frequency, large scale radio 
mapping for evidence of a plasma trail. 

The large scale, asymmetric r adio emission recently de- 
tected surrounding Cygnus X-3 feanchez- STu^etliil [2008b 
is also suggestive and we interpret it as dynamical interac- 
tion with the ISM. The overall morphology suggests that the 
source is moving in a roughly northern direction relative to 
the ISM. This should be followed up with VLBI proper mo- 
tion measurements of the XRB itself to test this hypothesis. 

4.2. Hot spots 

Given the estimate of the jet length l- ]ct from eq. ([8]l and the 
pressure inside the working surface, we can estimate the order 
of magnitude of the synchrotron emission from each hot spot 

as L v hs ~ e v —^p—. Once the plasma has moved out of the 
working surface of a jet, it expands (i.e., cools) adiabatically 
and reaches pressure equilibrium with the lobe. The pressure 
inside the hot spot is given by 



27r(aijot) c 
giving an approximate the hot spot luminosoity of 



^HS,nu- 



(34) 



(35) 



j 4 x 10 19 



ergs s 



Hz 



w 3/2 1/4 a l 



1/2 



/28r/? jet rc 



M 



V 25T+15 



i/4 



For typical microquasar parameters, the synchrotron cool- 
ing time rx for X-ray emitting electrons is longer than the 
adiabatic cooling time and the same expression should hold 
for radio as well as X-ray synchrotron emission. This flux 
is accessible with Chandra and in fact has bee n observed in 
a num ber of cases, most notably XTE J 1550 (ICorbel et alJ 
120021) . In the case of XTE J1550, the hot spot has been 
observed to move away from the XRB at mildly relativistic 
speeds, which fits in well with the picture developed here: If 
the jet power rises sharply, as indicated by the radio flares 
that precede the detection of the hot spots, the working sur- 
face will (a) brighten and (b) move outward, decelerating as 
it moves into undisturbed ISM. The distances where the hot 
spots are observed are consistent with where we would expect 
them from equation (0. 

4.3. Bow shocks and shells 

As has been dis cussed in the literature dGallo et al J 120051 : 
iRussell et al.ll2007h . and as is observed around radio galaxies 
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in clusters (e.g. [Birzan et al.l l2004) the source will sweep up 
a quasi-spherical shell during its initial inflation. This shell 
will emit thermal bremsstrahlung, potentially up to X-ray en- 
ergies, as well as line emission from the shocked plasma. 

For the early, quasi-spherical evolution and for stationary 
source, we refer the reader to the discussio n in the contex t 
of expandin g shells in galaxy cl usters (e.g. lHeinz et al.|[l998l) 
and binaries ( H einz et al.l2006aT) and concentrate on the emis- 
sion from the ISM swept up ahead of radio trails in the case 
of XRBs with appreciable space velocity. 

Given that LMXB space velocities top out at a few hun- 
dred km/s, we should expect temperatures up to about a keV, 
which should be detectable by X-ray telescopes (which is also 
the case in pulsar-wind nebulae). The radiative cooling time 
strongly depends on the ISM density, but assuming that the 
shock is fully radiative gives an upper limit on the total lumi- 
nosity from the swept up ISM of 



-^shock < nR, 



2 PISMWxrb 



(36) 



: 5 x 10 35 ergs s 1 • 



W 37 



P-11.5VH 1 /* 



spread over a region close to a degree in size (thus hard to 
detect with the typical field of view of current X-ray instru- 
ments). 

Longer wavelengths seem more suitable for finding these 
thermal sources (specifically, narrow band imaging of promi- 
nent shock-excited lines on degree scales, as weJJ_as broad- 
band imaging to look for free-free emission lRusselll l2008). 
Given that the estimated expansion velocity of the shell 
around Cyg X-l is of order a few dozen km/s, we would 
expect surface brightness and total power to be comparable 
with those detected in the Cyg X-l nebula. This suggests that 
asymmetric bow shock emission might be a suitable way to 



identify sources for low frequency follow up. 



5. CONCLUSIONS 

We proposed that jets from X-ray binaries (microquasars) 
leave behind trails of non-thermal, synchrotron emitting 
plasma as they move through the interstellar medium. 
LMXBs are more likely to leave such trails due to their longer 
life expectancy and, more importantly, due to the higher ex- 
pected space velocities, while HMXBs likely produce more 
stationary radio lobes. 

The total plasma volume deposited by XRBs over the his- 
tory of the Galaxy is comparable to the total disk volume and 
constitutes a non-negligible fraction of the halo volume. We 
argue that a fraction / m ; x of order a few percent of this plasma 
is mixed into the thermal ISM. 

The magnetic field thus released can easily provide the seed 
field for Galactic dynamos to produce the fields observed in 
spiral galaxies even under the most conservative assumptions. 
Since HMXBs turn on rapidly after an episode of star forma- 
tion, this mechanism should be important for galactic mag- 
netic field production and maintenance within a time frame of 
about 10 7 yrs from the first supernovae. 

We showed that radio emission from the plasma inside the 
trail should be visible primarily at low frequencies, since the 
radiative cooling time of the plasma at GHz frequencies lim- 
its emission to a roughly spherical region around the binary. 
LOFAR, MWA, LWA and other future low frequency, wide 
field instruments are ideally suited for searches of this emis- 
sion. LMXBs with high detected proper motions like XTE 
Jl 1 18+480 will be the best candidates for such a search. 

We would like to thank A. Merloni and C. Kaiser for help- 
ful discussions. S.H. acknowledges support from NSF astro- 
physics grant 0707682. 
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